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Inflation and Oscillations of Universe in 4D Dilatonic Gravity 
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We investigate the inflation of Universe in a model of four dimensional dilatonic gravity with 
a massive dilaton field <f>. The dilaton plays simultaneously the roles of an inflation field and a 
quintessence field. It yields a sequential hyper- inflation with a gr aceful exit to asymptotic de Sitter 
space-time, which is an attractor, and is approached as exp(—y/3A o6s ct/2). The time duration of 
the inflation is reciprocal to the the mass of the dilaton: A tinfi ~ m” 1 . The typical number of 
e-folds in the simplest model of this type is shown to be realistic without fine tuning. 
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I. INTRODUCTION 

We consider a specific 4D-dilatonic model of gravity 
(4D-DG) with a gravi-dilaton action 00 : 

Ag,* = ~YkS rf 4 aJ ^iI( M + 2A ° bSC/ ( $ ))> (!) 

where n = 8ttGn/c 2 is the Einstein constant and A obs 
is the observed value of the cosmological constant. In 
our model we have only one unknown function - the cos¬ 
mological potential t/($) - which has to be chosen to 
comply with all gravitational experiments and observa¬ 
tions in laboratory experiments, in star systems, in as¬ 
trophysics and cosmology. In addition it has to solve the 
inverse cosmological problem , namely, to determine [/(4>) 
that reproduces the time evolution of the scale parame¬ 
ter a(t) in accordance with the Robertson-Walker (RW) 
model of Universe 0 , 0 ). Further on we use cosmological 
units, in which the cosmological constant A obs equals one 
00 

The dilaton field 4) is a quantity, reciprocal to the di¬ 
mensionless gravitational factor g: <E> = 1/g. As usual in 
our scalar-tensor theory of gravity of Brans-Dicke type, 
instead of the Newton gravitational constant Gat, we de¬ 
fine the variable gravitational factor G e //( x) = Gn g[x) 
as a function of space-time coordinates x = {cf,x}. 

There exist different reasons for considering this model. 
For example, it appears in the five-dimensional Kaluza- 
Klein model, or, it can be produced in the low energy 
limit of string theory, by ignoring (yet unobserved) pos¬ 
sible higher dimensions. (See for details and for refer¬ 
ences 0.) The most physically relevant motivation for 
the use of a model with action 0 at present can be find 
in the attempts to develop a quantum version of Einstein 
gravity in four dimensional space-time |3|. 

Recently, a very promising discovery of the existence of 
a consistent quantum gravity in four dimensional space- 
time was made by Reuter 0. It confirms an early predic¬ 
tion by Weinberg that Einstein gravity in four dimensions 
may be ” asymptotically safe”, because of the existence of 
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a fixed point of the renormalization-group equations, sim¬ 
ilar to the one in 2 + e dimensions [g. Indeed, by making 
use of suitable Einstein-Hibert truncation, i.e., using ac¬ 
tion 0 with 4> = 1 and U (<E>) = 1 it has been proven that 
the four-dimensional possesses an ultra-violet fixed point 
and is asymptotically free. Then, as a result of the renor¬ 
malization, one obtains a running gravitational constant 
G and a cosmological constant A. In the case of quan¬ 
tum electrodynamics we have running electric charges, 
due to the cloud of virtual photons. In quantum gravity, 
a test body would be surrounded by a cloud of virtual 
gravitons. In contrast to electrodynamics though, the 
gravity is universal attractive and the effective mass of 
the body, as seen by a distant observer, is larger than it 
would be in absence of quantum effects. Hence, in quan¬ 
tum gravity we have an anti-screening effect, which en¬ 
tails the Newton constant and the cosmological constant 
becoming scale dependent quantities. As a result of the 
general relativistic invariance they become functions of 
space-time points. 

We prefer to describe the above phenomena using the 
dilaton field 4> in action 0 . This field encapsulates the 
space-time dependence of the corresponding parameters. 
The use of a cosmological potential C/(4>) of general form 
is known to be equivalent to the nonlinear gravity (NLG), 
described by a nonlinear (with respect to the space-time 
scalar curvature i?)Hilbert-Einstein Lagrangian L(R) in 
the action. (See 0 0 for more details and for refer¬ 
ences.) At present, the form of this Lagrangian is com¬ 
plete unknown. The renorm-group analysis of the non¬ 
linear terms in Reuter approach is still an open problem. 
It is clear that, in general, all coefficients in the different 
terms of the Taylor (or Laurent) series expansion of the 
function L(R) will become running constants, as a result 
of quantum effects. 

As shown in 00 , the use of the cosmological poten¬ 
tial U (4>), as an alternative description of these nonlinear 
terms, is extremely useful. Indeed, one is able to formu¬ 
late simple and clear physical requirements for this po¬ 
tential. For example, to have an unique physical de Sit¬ 
ter (classical) vacuum (dSV) in the theory, the function 
[/(4>) must have an unique positive minimum for positive 
values of the dilaton field 4>. To avoid the non-physical 
negative values of the dilation field $, the potential must 
increase to infinity for $ —► +0, etc. (See 0 .) 
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Here we use the simplest one-parameter potential: 

U{*) = & + W(< b) = $ 2 + A p -2 ($ _ I) (2) 

with such properties. The parameter is the dimension¬ 
less Compton length of the dilaton in cosmological units, 
defined by the equation p 2 = A obs l 2 , where 
is the usual Compton length of the dilaton, and is 
the mass of the dilaton. To comply with all known 
facts in gravitational physics, astrophysics and cosmol¬ 
ogy the dilaton mass must obey the constraint equation 
to 4 > 10 _3 eV. This yields an extremely small value of 
the dimensionless parameter < 10 30 110 

The cosmological potential 0 corresponds to a 
non-polynomial Lagrangian for the equivalent NLG: 
L(R\p tI> ) ~ —/(i?;p # ). It is described in parametric 
form by the equations: R = | p~ 2 (l/<b 3 — 4>) — 4$ and 
/ = | P~ 2 (3/<b 2 — 4> 2 — 2) — 2d) 2 , d> e (0,oo). One can 
see that the physical branch of the function f(R;p 4 ) has 
a Taylor series expansion 

f(R; P J=R + 2 + 1/3 (R + 4)(R + 8)pl+0(pl) 

and goes to the linear Einstein-Hilbert one f(R;p$) —> 
R + 2 (in cosmological units) in the limit p s —> 0, i.e., 
when —> oo. Hence, in this sense, our 4D-DG has the 
GR as a limiting case. 

As seen in Fig. Q for values of the scalar curvature 
\R\ <C 1 our model behaves like GR with standard New¬ 
ton constant Gjv ■ For R <C — 1 it behaves like GR with 
smaller effective gravitational constant, and for R 3> 1 - 
like GR with bigger effective gravitational constant: 

^ e //lji«-i ^ ^ e //l[Hi<ei ~ Gn < G e ff | R>1 - 


investigation of the effect of the presence of matter in 
Reuter approach to quantum gravity was started in JtJ|. 
Some general properties of 4D-DG model with matter 
were derived in [2j. In the present article we shell ignore 
the influence of matter for simplicity. Its presence yields 
a new important phenomena, which we intend to describe 
in a separate article. 


II. THE BASIC EQUATIONS FOR 
ROBERTSON-WALKER (RW) UNIVERSE IN 
THE 4D-DILATONIC GRAVITY 

Consider RW adiabatic homogeneous isotropic Uni¬ 
verse with 


ds 2 RW = dt 2 — a 2 (i)cMfc, 

where dl j? = + l 2 (dd 2 + sin 2 9)d<p 2 (for k = 

— 1,0,1). In our cosmological units the time t , and the 
scale parameter of Universe a(t) are dimensionless. 

For a free of matter Universe we obtain the following 
basic dynamical equations, governing the time evolution 
of the 4D-DG-RW model: 



= 


The dynamical equation for the dilaton <F is 


$ + 3-$ + U*($) =0 
a 


and 


( 3 ) 

( 4 ) 



FIG. 1: The form of the function /(R;p # ) in the equivalent 
nonlinear gravity for different values of the parameter p, = 
1/10, 1/20,1/30,'..., 0. 

In the general case, one must consider a 4D-DG model 
with a matter sector and a gravi-dilaton sector. The 


n*) = \v~ 2 ($ + | - 2^ (5) 

is the dilatonic potential jp). 

A. First Normal Form of the Dynamical Equations 

Introducing a new regularizing variable r by dt = Hdr 
and the notations X\ := H 2 > 0, X 2 := 4> > 0, X 3 := In a 
(where H = a/a is the Hubble parameter) we can rewrite 
the basic system © in standard normal form: 

i x=t(x) ' (6) 

where x,f £ are three-dimensional vector-columns 
with components {aq, X 2 , £ 3 } and {fi,f 2 l fs}, respec¬ 
tively, and 

fi = 4aq (Zx(x 2 ) - aq - ke~ 2x3 /2), 
h = x 2 (Z 2 (x 2 ) - xx - ke~ 2x3 ) , 

/3 = xi- 


( 7 ) 
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Here the quantities 


Zi($) = ^,*($) = Z 2 {$) + \vJQ, 
z 2 ($) = ![/($)/$ = +(8) 

are regarded as functions of x 2 = d 1 . 

Now one can derive another important relation of a 
pseudo-energetic type by using Eq. 0 . The qualitative 
dynamics of its solutions is determined by the function 

v = U 2 + v($). (9) 


This function plays the role of a (non conserved) energy¬ 
like function for the dilaton field < f > in the 4D-DG-RW 
Universe and obeys the equation 



( 10 ) 


As seen from this equation, the quantity 77 defines a 
Lyapunov function in a generalized sense in our model of 
Universe, filled only with vacuum energy. In this case the 
parameter 77 ( 7 ") is a monotonically decreasing function of 
the time r. One can show, that the same relations are 
valid for the case of ultra relativistic matter, when the 
trace of energy-momentum tensor equals zero and the 
function 77 is a Lyapunov function for the corresponding 
generalization of the normal system of ordinary differen¬ 
tial equations © 13- 

Ail solutions to the system 0 must cross the surfaces 
of level 77 = const of this function and go into its inte¬ 
rior. The form of this two-sheeted surfaces is described 
explicitly by the equation: 



FIG. 2: Level surface of the Lyapunov function rj - Eq. 0 
for the potential V($) - Eq. @ and k = +1. 
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H = dy/ 77 -V"+$A ± /( 77 -V) (77 —U+2$A). ( 11 ) 

Here A = -|U — k-%?. A typical geometrical form of a 
level-surface of the function 77 in absence of matter for 
the simple potential V = V($) 0 is shown in Fig. 0 - 
for k = 1, and in Fig. 0 for k = —1. For k = 0 the 
function tfTTj) does not depend on the scale factor a and 
the corresponding surface is a simple cylinder. 

One can see two main specific features of the behavior 
of this surface: 

1) At small values of the RW scale factor a ~ 0 the form 
of the Lyapunov surface depends strongly on the sign of 
3-space scalar curvature: For k = +1 we have a semi¬ 
infinite curved cylinder with a closed bottom, which goes 
down to values H = 0 in ’’vertical” direction. In contrast, 
for k = — 1 we have on both sides an infinite curved 
cylinder which goes up to values H = 00 in ’’vertical” 
direction. At the end, for k = 0 we have a simple straight 
’’horizontal” cylinder. 

2) For all values k = 0, ±1 and a —> 00 we have in¬ 
finitely long “horizontal” tube along the a-axes. The so¬ 
lutions that enter the inner part of the {a,Q,H}- space 


FIG. 3: Level surface of the Lyapunov function 77 - Eq. 0 
for the potential U(<f>) - Eq. 0 and k = — 1. 

through the Lyapunov surface, in this area go to de Sitter 
asymptotic regime winding around the axis of the tube 
and approaching it for t —> 00 . 

The analytical behavior of the solutions in the com¬ 
pletely different regimes at a ~ 0, and at a 1, as well 
as the transition of given solution to asymptotically de 
Sitter regime will be described below in more detail. 

B. Second Normal Form of the Dynamical 
Equations 

We need one more normal form of the dynamical 
equations to study the behavior of their solutions for 
X\,X 2 —7 00 . The standard techniques for investigation of 
the solutions in this limit show that the infinite point is a 
complex singular point, and one has to use the so called 
cr-process |8| to split this singular point into elementary 
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ones. This dictates the following change of variables: 


° —2 —4 —4 -1 . ^ /-> 

Xl = 16 P * * 9 ' X2 = 9 1 t= ^ P ' s • &, 


( 12 ) 


III. GENERAL PROPERTIES OF THE 
SOLUTIONS IN THE 4D-DG RW UNIVERSE 


1. Properties of the Solutions in a Vicinity of dSV 


which transforms the equations © with right hand sides 
© into a new system: 

g' = gV, z' = zZ, <d’ = z 2 g 2 (13) 


where the prime denotes differentiation with respect to 
the variable X 3 and 


V(g,z,a;pl,k) 


Z(g,z,a;pl,k) 


1 - ^z 4 g 5 w(g) 


%-) g 3 z 4 ( 1 —3.g 


V 0 (g,z)+O 2 (pi) 

\z i g 7 v, g {g) 

2 


3 ' 4 p *' 9 V|- 


2 V P) g ' z 'C 
Z 0 (g,z)+O 2 (pl) ■ 


(14) 


Here we have introduced the functions v(g) = -yP^U(i) 
and w(g) = -yp 2 Wig) which do not depend on the param¬ 
eter and are simply related: w = 3/2 g ~ 2 J® g 3 dv. The 
functions 2? and 2 are defined according to the formulae 


= ri(lng) = d( In z) 

d(lno) ’ d(lna) ’ 


The representation da shows that one can develop a 
simple perturbation theory for the highly nonlinear sys¬ 
tem da using the extremely small parameter p 2 as a 
perturbation parameter. 

From the dynamical equations ®, © and m , one 
easily obtains the following contour-integral representa¬ 
tion for the number of e-folds J\f and for the elapsed time 

AtAr = 73 P * 0Ar: 


Tfin 

N(pl,k)= J H 2 {r)dT= J 


dg/g 


v{g,z,a-,p\\,k)' 


(16) 


and 




z 2 gdg 


'Dig, z , a-,p 2 , k) ’ 


(17) 


Here, a start from some initial (in) state of the 4D-DG- 
RW Universe, followed by a motion along the contour 
C,£" (determined by corresponding solution of system 
Ill3l) h and an end at some final (fin) state, are assumed. 


Let us first consider the simplest case when k = 0. 
The system © in this case splits into a single equation 
for X 3 , which is solved by the monotonic function X 3 (r) = 
£3 + f xi(r)dT, and the independent of X 3 system 

4~xi = 4xi (Zi(x 2 ) - £ 1 ), 
ar 

4-X 2 = X 2 (Z 2 (x 2 ) - Xl) . (18) 

Now it is clear that the curves X\ = Zi{x 2 ) and 
x\ = Z 2 {x 2 ) are the zero-isoclinic lines for the solutions 
of (1181) . These curves describe the points of local extrema 
of the functions X\ (r) and x 2 (r), respectively, in the do¬ 
main X\ 2 > 0. Because of the condition U(<h) > 0 and 
the existence of a unique minimum of the cosmological 
potential, these lines have a unique intersection point - a 
de Sitter vacuum (dSV) state |2‘] with X\ = 1/3, x 2 = 1. 
This singular point represents the standard de Sitter so¬ 
lution, which in usual variables reads 

H = l/-v/3, $ = 1, a(t) = ao exp(t/v / 3). (19) 

One can see the typical behavior of the solutions of 
the system m in the domain 2:1,2 > 0 , together with 
the curves Z[l] := {x\ = Z 1 ( 0 : 2 )}, and Z{ 2] := (ah = 
Z 2 (x 2 )}, in FigQJ where the corresponding phase portrait 
is shown for the case of the potentials © and p^ = 1/4. 



FIG. 4: A typical phase portrait of the system ifflfll . The 
parts of solutions with the same color are covered by the Uni¬ 
verse for equal r-intervals. 

Consider the solutions of the system m of the form 
xi t 2 {r) = Xi ]2 + Sx 1 i 2 (t) that are close to dSV, i.e., with 
|<5ari,2 (t)I 1. Using the relations 

V(1)=0, V,*(1) = 0, U, $$ (l)=p; 2 ; 

U(l) = l, [/,*(!)= 2, C7,**(l) = 2^1 + |p" 2 V(20) 
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which describe the general properties of de Sitter vacuum 
state (with <f> = 1) in 4D-DG Q, one obtains in linear 
approximation 


/ ^0 \ 

5x\(t) = 5x\l exp —— t ] cos(w # f), 


3 pjx\ / \/3\ , * , ^ 

Sx 2 {t) = - exp — —t cos 


2 <Pi +! 


V3,' 


Sx- 3 (t) =Sx° + \/3p^Sxi exp (—— f] x 


(cos (ui^t-lpa) - COS (ipa)) , (21) 


where 


= \Jp~i 2 - 3 / 4 > 

, 2^3 , 2y/3 

tan ip^ = —— , tanyv =-— (22) 

5 o 

Sx ° = fei(O) is the small initial amplitude of 5xi(t), and 
the solution for the deviation Sx 3 (t) has been added for 
completeness and for later use. 

The frequency is a real and positive number if p s < 
2/V3- According to the estimate p # < 10 _3 °, in 4D- 
DG we have > 10 3 ° in cosmological units, or <x> # > 
10 3 GHz in usual units, and ~ 7r/2 ss — ip a . 

Thus, we see that the de Sitter vacuum is a stable 
focus in the phase portrait of the system 11811 . For k = 0 
all solutions of this system that lie in a small enough 
vicinity of dSV oscillate with an ultra-high frequency 
<!221l . and approach dSV in the limit t —> oo. 

One can generalize this statement for the case of ar¬ 
bitrary k = 0, ±1 H]. This turns out to be true even 
in the presence of nonzero matter energy-density terms. 
Then we see that one important general prediction for 
the 4D-DG-RW Universe is the existence of ultra-high 
dilatonic oscillations with frequency (PI) in 3-spaces 
with any curvature and in the presence of any kind of 
normal matter 0. 

If P* > 2 /a/ 3, i.e., if m t ~ 10 _33 eV (as in inflation 
models with a slow-rolling scalar held and in quintessence 
models 0,0,0), the above ultra-high dilatonic oscil¬ 
lations do not exist in the 4D-DG-RW Universe. In this 
case, the frequency u>$ becomes imaginary and, instead 
of a stable focus, we have an unstable saddle point in the 
phase portrait of the system ©■ Such a situation was 
considered first in 0 in a different model of nonlinear 
gravity based on a quadratic with respect to scalar cur¬ 
vature R gravitational Lagrangian, but with some addi¬ 
tional terms that originate from quantum fluctuations in 
curved space-time 0. These additional terms are related 
to the Weyl conformal curvature and therefore vanish in 
the case of RW metric, but yield an essentially differ¬ 
ent theory in other cases. Therefore, for RW Universe 
the model described in [u| is equivalent to 4D-DG with 
a quadratic in the held >!> cosmological potential. Such 


a potential is non-physical, because it admits negative 
values of the dilaton held H>, i.e., of the gravitational 
constant. 

An immediate consequence of the above consideration 
is the existence of ultra-high frequency oscillations of the 
effective gravitational factor G e // = G N /$, accompa¬ 
nied with an extremely slow exponential decrease of its 
amplitude ~ exp lV3A obs ct^j (in usual units). One 
obtains H 2 /Hq = Ua and SHq = 1 • Hence, at the 

present epoch with ss we have 8 x\ « g. Then the 
second equation of the system CD gives 


/q / /o \ 

g(t) « 1-P* — exp( ——t J cos , (23) 


where we have introduced a dimensionless gravitational 
factor, g(t ) = G eff (t)/G N = 1 /$(£)■ 

Because of the extremely small amplitude p # < 1(U 30 , 
these variations are beyond the possibilities of present- 
day experimental techniques. 

In contrast, the oscillations of the Hubble parameter 
H have a relatively big amplitude SH 0 = \fbx o ~ 0.4, 
and the same huge frequency as the oscillations of 
gravitational factor. It is very interesting to find possible 
observational consequences of such a phenomenon. 

High frequency oscillations of the effective gravita¬ 
tional factor were considered first in the context of Brans- 
Dicke field with a BD parameter u> > 1 in 0. These os¬ 
cillations were induced by an independent inflation field, 
but the analysis of the existing astrophysical and cos¬ 
mological limits on the oscillations of G e //(t) is appli¬ 
cable to our 4D-DG model as well. The conclusion in 
{0 is that the oscillations in the considered frequency- 
amplitude range, being proportional to g/{gH ), do not 
affect the Earth-surface laboratory measurements, Solar 
System gravitational experiments, stellar evolution, and 
nucleosynthesis, but can produce significant cosmological 
effects because the frequency is very large and the Hubble 
parameter is small (in usual units). It can be seen explic¬ 
itly from Eq. <ED that this is precisely what happens in 
4D-DG, although in it the oscillations are self-induced. 

As stressed in 0 , despite the fact that the variations 
of the type ED have extremely small amplitudes, they 
can produce significant cosmological effects because of 
the nonlinear character of gravity. In the case of free 
of matter Universe the 4D-DG version of the formula, 
analogous to the one in the first of references 0, is 


h= 1 - 9 -± 

2 9 




:<JU 



(24) 


Being a direct consequence of Eq. 0 , this formula 
shows that, after averaging of the oscillations, the term 
g/g has a non-vanishing contribution because it enters 
the Hubble parameter ED in a nonlinear manner. A 
more detailed mathematical treatment of this new phe¬ 
nomenon in 4D-DG is needed to derive reliable conclu¬ 
sions. The standard averaging techniques for differential 
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equations with fast oscillating solutions and slowly devel¬ 
oping modes seem to be the most natural mathematical 
method for this purpose, but the applications of these 
techniques to 4D-DG lie beyond the scope of the present 
article. 


IV. INFLATION IN 4D-DG-RW UNIVERSE 

Having in mind that: 1) the essence of the inflation 
is a fast and huge re-scaling of the Universe @01 , and 
2 ) the dilaton is the scalar field determining the scales 
in Universe, it seems natural to relate these two funda¬ 
mental physical notions instead of inventing some specific 
“inflation field”. In this section we show that our 4D-DG 
model indeed offers such a possibility. 


A. The Phase-Space Domain of Inflation 

As seen from the phase portraits in Fig. for val¬ 
ues H > H cr n and > $ CT .j t , the ultra-high-frequency- 
oscillations do not exist. The evolution of the Universe 
in this domain of the phase space of the system © re¬ 
duces to some kind of monotonic expansion, according 
to the equation 4-x 3 = H 2 > 0. We call this expansion 
an inflation. As we shall see, it indeed has all needed 
properties to be considered as an inflation phenomenon. 

The transition from inflation to high-frequency oscil¬ 
lations is a nonlinear phenomenon, and we will describe 
it in the present article very briefly. Here our goal is to 
have some approximate criteria for determining the end 
of the inflation. It is needed for evaluation of the basic 
quantities that describe the inflation. 

As seen from Eq. m, the amplitude, <5$o, of the 
oscillations of $ is extremely small compared with the 
amplitude SH 2 of the oscillations of H 2 : <5d>o < 

IpJ^/pI + | 8Hq. An obvious estimate for the ampli¬ 
tude 5Hq is SHq < 1/3 (= H 2 ). Then, for < 2/\/3 
(which is the condition for existence of oscillations), we 
obtain <5$o < SHq < 1/3. The last estimate is indeed 
very crude for the physical model at hand, in which 
< 1(T 30 . This consideration produces the constraint 
Hcrit -\/2/3 and § C rit <4/3, but, taking into account 
the extremely small value of p # , we will use for simplicity 
the very crude estimate H crit , <1> crit ~ 1. 

Now it becomes clear that the study of the inflation 
requires considering big values of the variables Xi <2 , he., 
using the second normal form m of the dynamical equa¬ 
tions. 


B. The Behavior of the Solutions Near the 
Beginning 


1. The Case k = 0 

Let us consider first the case k = 0. From Eq. m one 
obtains the simple first order equation 


(h_z 5 g 6 v tg {g) 
dg 


4 v(g,z;pZ) 


(25) 


with 


V(g,Z',Pj = 1 - 3 z 9 u{g) = 


h-g- g Ag)- 


4 Pt 


<i> \ 3 4 


g 3 Z * = V 0 (g,z)+O 2 (pl), (26) 


where u{g) = *fp\U{l/g). 

One can easily prove that the solutions of Eq. m have 
important general properties for small values of g. For 
the potentials © we obtain 


(g,z) = -^-g 5 (l- g 2 ) 


(27) 


and 


T> 0 (g,z)=i- : jg 3 (1 -g 2 ) 2 ■ 


(28) 


Hence, lim g ^o T > 0 (g, z) = 1. Then, taking into account 
the leading terms, we obtain from Eq. tm U3, U3- 
(1251) . and cu the following results in the limit g —* 0 : 

a(g) ~ g, 


z(g)=z 0 - —4g 5 + °(g 6 )> 
t(g)= 2 J±z 2 g 2 + 0(g 3 ). 


(29) 


When solved with respect to the time t, these formulae 
give 

1/2 

I +0(t), 


a(t) ~ — 

\P* 


t 


1/2 


Zo \ 2 , 


z(t) = z 0 - 


15 


Al/ 

■ 2 ''* / 


5/2 


0{t 3 ). (30) 


2. The Case k = ± 1 

Now, the equations do not split into a two- 
dimensional subsystem and an independent differen¬ 
tial equation and must be considered as a true three- 
dimensional system. The expansion of its solutions in a 
Taylor series with respect to the variable a is 

g{a) = gi a + 0(a 2 ), 

z{a) = z 0 - ^ kpl gf a 2 + 0(a 3 ), 

t ( a ) = f^ p * z o 9i a 2 + 0(a 3 ), 


(31) 
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with an arbitrary constants g± and zo- When solved with 
respect to the time t, these formulae give 

a{t) = z^gf 1 

g(t) = zo 1 + °( t ) 

z(t) = zo-k(^j t + 0(t 3 ' 2 ). (32) 

We see, that the value of the sign k = 0, ±1 of the 3- 
space curvature influences only the form of the functions 
z(t) and z(a ) for small values of the time t, and of the 
scale parameter a, respectively. These differences in the 
behavior of the functions z(t) and z(a ) for different val¬ 
ues of the parameter k are consistent with the pictures, 
shown in Fig. 0 and Fig. [3 


3. Some General Conclusions 

The above formulae <!291l . cp, on . and on show for 
any value of parameter k = 0, ±1: 

1) the existence of Beginning, i.e., the existence of a 
time instant t = 0 at which a(0) = 0; 

2) the asymptotic freedom of gravity, i.e., the zero 
value of gravity at the Beginning: g(0) = 0; 

3) the finiteness of the time interval needed for reaching 
nonzero values of a and g , starting from the Beginning; 

4) the high precision constancy of zft) at the Beginning 
for k = 0 and the dependence of linear term in the Taylor 
series expansion of the function z(t) for small values of 
cosmic time t on the sign k = ±1 of the 3-space curvature. 

5) the behavior of the RW scale factor aft) for small 
t , similar to its behavior in GR in the presence of radia¬ 
tion. Hence, one can conclude that at the Beginning the 
dilaton plays a role, similar to the role of a radiation. 

The asymptotic freedom of gravity at the Beginning in 
4D-DG-RW Universe is consistent with the qualitative 
results in 0], It leads to some sort of an initial power- 
law expansion. The number of e-folds A fft) —> — oo as 
t —> +0 like lnt, since Afft) = lna(f). 

The general conclusion is that the potentials ©, © 
do not help to overcome the initial singularity problem in 
the free of matter 4D-DG-RW Universe with k = 0, ±1. 
However, we would like to emphasize that: 

1) 4D-DG with the simplest potentials ©, © may 
not be applicable for times smaller then the Planck time 
tpi ~ 10~ 44 sec, because it may be a low-energy the¬ 
ory, which ignores some essential quantum corrections 
in gravi-dilaton sector. (See 3 for derivation of 4D- 
DG as a low-energy limit of the string model.) But one 
may expect 4D-DG to be valid after some initial time 
instant, t in ~ t# = h/fm^c 2 ). If m# <C M P i, we will 
have ti n tpi and our results for the case under consid¬ 
eration may have a physical meaning, leaving open the 
initial singularity problem. 


2) Another possibility to overcome the initial singular¬ 
ity problem is to introduce some kind of matter in the 
Universe. The asymptotic freedom of gravity at the Be¬ 
ginning shows that the gravitational attraction goes to 
zero for small values of the scale parameter a. Then, the 
existence of some amount of usual matter will lead to 
a positive pressure in the Universe and to correspond¬ 
ing repealing forces, which will not be compensated by 
gravity at small values of the scale factor a. The prelim¬ 
inary numerical investigation showed that for a 4D-DG- 
RW Universe with normal matter (dust, perfect fluid, 
and/or radiation) one indeed has a bouncing solutions 
aft) ^ 0 without initial singularities. 


C. The Inflation and the Number of E-Folds 


Our general analysis of the solutions for 3D-DG RW 
model of Universe without matter shows, both analyti¬ 
cally and numerically, that its qualitative behavior does 
not depend essentially on the sign of the 3-space cur¬ 
vature. Therefore, for simplicity, we consider in what 
follows only the case k = 0. 

According to Eq. ITUI l and Eq. HJ one can represent 
the RW scale factor aft) in the form: 


aft) = gft)exp(AN'ft)), 
where the re-normalized number of e-folds: 

AAf f dg 1-V{g,z,a;pl,k) 

J 9 ^(g,z,a;pl,k) 

c/‘ n 


(33) 


(34) 


is a finite quantity in the entire time interval t £ [0,oo). 
Obviously, AjV(fW) = Afft^) at the special time in¬ 
stants iW, i = 0,1,..., when gft) reaches its dSV value, 
i.e., when gft W) = 1 , for example, in the limit t —> oo. 

The phase portrait of Eq. m and the time- 
dependence of the dimensionless gravitational factor gft) 
for the potentials ©, © are shown in Fig. [U] and Fig. [U] 
respectively. From Fig. [U] we see that one can define 
analytically the time of duration of the initial inflation 
as the time spent by Universe from the Begin¬ 
ning to the first time instant t^°\ when the gravitational 
factor g(t^) = 1. In addition, we see in Fig. |© that this 
time interval is finite and has different values for different 
solutions of Eq. ©2). 

In Fig. © and Fig. one can see a new specific fea¬ 
ture of 4D-DG: the solutions may enter many times the 
phase-space domain of inflation and the function gft) os¬ 
cillates around its dSV value g = 1 with a variable period. 
Between two successive maxima of gft), the squared log¬ 
arithmic derivative H 2 of aft) has its own maxima, just 
at the already defined time instants In the vicinity 
of each maximum of H 2 , the function a(t) increases very 
fast - like exp(const x fl*), with const > 0, and parameter 
H > 2. Therefore, we call such inflation, which is much 
faster than the usual exponential de Sitter inflation, a 
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FIG. 5: The phase portrait of Eq. 1251 . The thin line shows 
the zero line Z[2] of the denominator V. 



FIG. 6: The dependence of dimensionless gravitational fac¬ 
tor g on cosmic time t for solutions of Eq. (1251 . 

hyper-inflation. Hence, in 4D-DG, we have some sort of 
successive hyper-inflations in the Universe. For simplic¬ 
ity, we define the time-duration A t^f t = t W — 

i = 0,1,... (by definition A= t®) of each of these 
periods of hyper-inflation as the time period between two 
successive dSV values of the function g(t ), although the 
hyper-inflation itself takes place only around the max¬ 
ima of the function H‘ 2 (t). The corresponding number 
of e-folds is N$ f v It is clear that A f^ t is a decreasing 
function of the number i. The inflation can be considered 
as cosmologically significant only if for some short total 
time period At‘°^ a / = A t^fi, the total number of 

e-folds A QZjf = (T Aexceeds some large enough 

number AC. It is known that one needs to have AC > 60 to 
be able to explain the special-flatness problem, the hori¬ 
zon problem, and the large-scales-smoothness problem in 
cosmology j9]. 

Fig- 0 illustrates both the inflation and the asymptotic 
behavior of the function A f(t) for t —> oo. The small 
oscillations of this function are “averaged” by the crude 



FIG. 7: The dependence of the number of e-folds Af on the 
cosmic time t for solutions of Eq. ifHt . The straight line gives 
an example of an asymptotic line for some solution. 

graphical abilities of the drawing device, and for large 
values of time t we actually can see in Fig. [3 only the 
limiting de Sitter regime m, when the averaged func¬ 
tion (A f{t)) = (AAf(t)}. In this regime, for t —> oo, we 
have an obvious asymptote of the form (A f(t)) ~ A T + Ht 
for the averaged with respect to the dilatonic oscillation 
m function (Af(t)). We accept the constant AC as our 
final definition of total number of e-folds during inflation: 

A 7 = lim ((AC(i)) — Ht) = lim ((AAC(i)) — fit ) . (35) 

It is clear that this integral characteristic describes pre¬ 
cisely the number of e-folds due to the true inflation in 
the 4D-DG-RW model, i.e., during the fast initial ex¬ 
pansion of the Universe, as a new physical phenomenon. 
To obtain this quantity, one obviously must subtract the 
asymptotic de Sitter expansion from the total function 
(AC(t)). It can be shown that one has to apply the same 
definition to the solutions in the general case of arbitrary 
values of the space-curvature parameter k and physically 
admissible matter energy densities e(a), since they have 
the same asymptotic behavior. Hence, in all cases the 
solutions with a cosmologically significant inflation must 
have values J\f > 60. 

As we see in Fig. [3 the total number of e-folds AC 
decreases, starting from the solution “4” (with = 0.2) 
to the solution “3” (with z o = 0.7) and to the solution 
“2” (with zo = 2.7), reaches a minimum (sa 1) for some 
initial value Zq, and then increases for the solution “1” 
(with zq = 7), and for solutions with larger values of Zo- 
Thus, we see that inflation is a typical behavior for all 
solutions of Eq. (ESI), and that most of them have large 
values of Af. (In Fig. [3 we show only solutions with small 
values of AC that are close to its minima.) 

The value of the constant Af depends on the closeness 
of the given solution to the zero curve Z[2] of the denom¬ 
inator of the integrand in the right-hand side of ifTTfli . 
The equation V = 0 can be explicitly solved. Its solu¬ 
tion reads z = g 5 u(g )) . The corresponding curve 

is shown both in Fig. Q] and Fig. [3 
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The value of the quantity A f(t) increases essentially 
each time when the solution approaches this curve. The 
corresponding increment of A f(t) remains finite, even if 
the solution crosses this curve, although in this case the 
denominator in the integrals El), or reaches a zero 
value. 

Indeed, taking into account that the points where 
the solution crosses the curve Z[ 2] are extreme points 
of the function g(z), and using the expansion Az = 
1 /2 

(—+ Oi(Ag) in a vicinity of such a point, 
we easily obtain 


AAf = z- A g- 6 u~\ -—Ag + O^Ag), 

V 9V,g 

At= —=z~ 2 g~ 4 u~ 1 J——Ag + 0 1 (Ag) (36) 

V O y gv,g 

for potentials u(g) and v(g) of the most general type. 
Here the values of the coefficients are taken on the curve 
V = 0. 

Combined with the previous results, this proves that 
the time intervals of inflation are finite for all values 
of i = 0 , 1 ,.... 

Thus we see that an essential increase in A/(f) is ac¬ 
cumulated when the solution becomes close and parallel 
to the curve Z[2\. This is possible both for big values of 
g and small values of 2 (the solution “4”), or for small 
values of g and big values of z (the solution “1”). All this 
results in a big value of J\f. Only a small fraction of the 
solutions (like the solutions “2” and “3”) stay all the time 
relatively far from the line T> = 0, and therefore acquire 
relatively small number of e-folds, 1 < A/" < 2 — 3. The 



FIG. 8: The dependence of the number of e-folds A f on the 
initial value zo of the solutions of Eq. (1251 . 


the following expression, obtained by a numerical analy¬ 
sis of the solutions of Eq. m- 

A T(z 0 ) « 2.552 z- 3/10 + 0.090 z® /3 . (37) 

Another important numerical observation is the depen¬ 
dence of the relation A/"(zo,p # ) on the value of the param¬ 
eter p 9 . As seen in Fig. El there is a good convergence of 
the corresponding family of curves, when —> 0. Having 
in mind that in the realistic case < 10” 30 , we see that 
the inflation in 4D-DG is a robust property with respect 
to variations of the mass of the dilaton. The position of 
the minimum of the number of e-folds is at zq ~ 1.6 and 
seems to be independent on this mass. The value of the 
minimum is ^ 2.5 in the limiting case p 4 —> 0. 



FIG. 9: The dependence of the number of e-folds A f on the 
initial value zo of the solutions of Eq. G3. The curves for dif¬ 
ferent values Pq = 1/50,1/100,1/150,1/200, and 1/250 are 
shown. As seen in FigED for Pq, —> 0 the family of curves 
converges to a limiting curve, which bounds this family from 
above. 

This qualitative consideration, combined with the 
structure of the phase portrait shown in Figures 0 and EJ 
not only explains the universal character of the inflation 
in 4D-DG, but gives us a better understanding of its ba¬ 
sic characteristics. In particular, we see that we do not 
need fine tuning of the model to describe the inflation as 
a typical physical phenomenon. 

Using Eqs. OB. (TT7I) . and to . we obtain 

4n/i(p«) = ^=p.e£,/i( o) + o(p 3 j, 

where 


dependence of the number of e-folds on the initial value 
Zq is shown in Fig. [S] Here one can see that the most of 
the solutions of Eq. m have big values of J\f. One can 
reach Af ~ 50 — 100 for zo ~ 10—15, i.e., without any fine 
tuning. An approximate form of this dependence (at a 
few percent level of accuracy) for = 1/250 is given by 



dg/g 

v(g,zioy 


and 


0 


0) 
inf l 


( 0 ) 


/ 




z 2 gdg 

V(g,z;0y 


(38) 
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are independent of the parameter p Taking into ac¬ 
count the extremely small physical value of this param¬ 
eter, one can conclude that the higher order terms in 
^"infi (p$) an d A are not essential. Neglecting 
them, we actually ignore the contribution of the term 
<h/3 in the functions Z\p, (0, or the corresponding term 
<f> 2 , in the cosmological potential 17(d)) ©. It is natu¬ 
ral to ignore these terms in the domain of inflation, be¬ 
cause they are essential only in a small vicinity of dSV. 
In the function W(d>), we have a term that dominates 
for d> — 1 p 2 , having a huge coefficient ~ p~ 2 . Phys¬ 
ically this approximation means that we are neglecting 
the small pure cosmological constant term in the cosmo¬ 
logical potential and preserve only the terms, which are 
proportional to the mass of the dilaton. 

The relation AtV fl (j%) ~ ^0^,, written in phys¬ 
ical units, reads 

E»At%~h- l0#(O). (39) 

It resembles some kind of a quantum “uncertainty rela¬ 
tion” for the rest energy E # of the dilaton and the time of 
inflation and probably indicates the quantum character 
of the inflation as a physical phenomenon. 

More important for us is the fact that Eq. shows 
the relationship between the mass of the dilaton, m,, and 
the time duration of the inflation. Having large enough 
mass of the dilaton, we will have small time duration of 
the inflation. This recovers the real meaning of the mass 
to # as a physical parameter in 4D-DG, and gives possi¬ 
bility to determine it from astrophysical observations as 
a basic cosmological parameter. 

V. CONCLUDING REMARKS 

In this section, we discuss some open problems in 4D- 
DG. 

The main open physical problem at the moment seems 
to be the precise determination of the dilaton mass. The 
known restriction on it is too weak. It is convenient 
to have a dilaton d> with mass m 4 in the range 10” 3 - 
1CT 1 eV. In this case, the dilaton will not be able to de¬ 
cay into the other particles of the Standard Model, since 
they would have greater masses 0- On the other hand, 
in 4D-DG we do not need such a suppressing mechanism 
since a direct interaction of the dilaton 4> with matter of 
any kind is forbidden by the weak equivalence principle 
i. This gives us the freedom to enlarge significantly the 
mass of the dilaton without contradiction with the known 
physical experiments. One of the important conclusions 
of the present article is that is related to the time 
duration of the inflation. One is tempted to try a new 
speculation - to investigate a 4D-DG with dilaton mass, 
to 4 , in the domain 100 GeV-1 TeV, and dilaton Comp¬ 
ton wave length, Z 4> , between 10~ 18 and 10~ 16 cm. In 
this case the time-duration of the inflation, , will be of 


order of 10~ 28 sec; the dimensionless dilaton parameter, 
Pj,, will be about 10~ 45 , and the ultra-high frequency, 
of dilatonic oscillations during de Sitter asymptotic 
regime will be approximately 10 19 GHz. Such new val¬ 
ues of the basic dilaton parameters are very far from the 
Planck scales. They seem to be accessible for the parti¬ 
cle accelerators in the near future, and raise new physical 
problems. 

Another open problem is the justification of quantum 
corrections and the exact form of cosmological potential. 
More general potentials than © where introduced in @], 
and further justification and verification is needed. One 
must take into account the possible additional corrections 
in order to obtain a correct physical description in the 
early Universe. 

Important open problems in the development of a gen¬ 
eral theoretical framework of 4D-DG are the detailed 
theory of cosmological perturbations, structure forma¬ 
tion, and possible consequences of our model for the 
CMB parameters. The properties of the solutions of the 
basic equation for linear perturbations 5 4) differ essen¬ 
tially from the ones in other cosmological models with 
one scalar field. This equation yields a strong “cluster¬ 
ization” of the dilaton *1) at very small distances 0. 
Actually, the equation for dilaton perturbations shows 
the existence of ultra-high frequency oscillations, de¬ 
scribed above, and non-stationary gravi-dilaton waves 
with length < 10~ 2 cm. Such new phenomena can¬ 
not be viewed as a clusterization at astrophysical scales. 
Thus, their investigation as unusual cosmological pertur¬ 
bations is an independent interesting issue. For example, 
it is interesting to know whether it is possible to consider 
these space-time oscillations as a kind of dark matter in 
the Universe. 

A more profound description of the inflation in 4D-DG, 
both in the absence and in the presence of matter and 
space-curvature, is needed. It requires correct averaging 
of dilatonic oscillations. 

We intend to present the corresponding results else¬ 
where. 
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